Competing orders in FeAs layers 
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Using the unrestricted Hartree-Fock approximation and Landau theory we identify possible phases 
competing with superconductivity in FeAs layers. We find that close to half-filling the transition from 
the paramagnet to the magnetically ordered phase is first-order making anharmonicities relevant and 
leading to a rich phase diagram. Between the already known one dimensionally modulated magnetic 
stripe phase and the paramagnet we find a new phase which has the same structure factor as the 
former but in which magnetic moments at nearest-neighbor sites are at right angles making electrons 
to acquire a non trivial phase when circulating a plaquette at strong coupling. Another competing 
phase has magnetic and charge order and may be stabilized by charged impurities. 
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Exotic superconductivity often appears in a region 
where a paramagnetic (PM) phase competes with ordered 
phases. Examples are the heavy fermion compounds[l(, 
organics[2j, Ijs, Jj], vanadium bronzes [5[ and barium 
bismuthates[6|, |7|. A competing phase is believed to 
play an important role in Cu based high temperature 
superconductors a |3, [lfj, [ill |l2|, |l3|, 14 1 probably with 
charge[15lll6lll7| and magnetic order 15 1 although there 
are other more exotic proposals [l2| [l3|, \1M- Alterna- 
tively one can take those as secondary effects and con- 
sider the su per conductor as emerging from the Mott 
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insulator [18|, [19|, |20|, |2l| or from the Fermi liquid with 
antifcrromagnetic correlations [22j . One important ques- 
tion is if the competing phase is only detrimental for 
superconductivity[23[ or if the vicinity to an ordering in- 
stability helps superconductivity [9. Ho. l2fj. |22{| . 

In the recently discovered Fe-based superconductors, 
the Coulomb interaction is believed to be relatively 
weaker than in cuprates making a mean-field analysis 
more likely to lead to a quick identification of the pos- 
sible competing phases. Indeed first principle computa- 
tion have predicted a magnetic stripe (MS) phase with 
ordering wave vector (w, 0) (we define wave vectors for an 
isolated FeAs layer with the primitive vectors connecting 
nearest Fe sites and Fe-Fe distance a = 1) which later 
has been found by magnetic neutron scattering 24j . 

Here we extend the analysis using the unrestricted 
Hartree-Fock approximation in a model Hamiltonian 
which keeps only two orbitals supplemented with a Lan- 
dau theory analysis. This allows more flexibility than 
first principle computations to identify possible compet- 
ing phases at the expense of more (yet) unknown parame- 
ters. Surprisingly, we find a new magnetic phase in which 
magnetic moments at nearest neighbor sites are at right 
angles which we term orthomagnetic (OM). This is inter- 
esting in its own because well formed magnetic moments 
usually lead to phases in which magnetic moments are 
either mutually parallel or antiparallel[25(. We identify 



another competing phase which has both charge and spin 
order reminiscent to the stripe phases of cuprates. This 
phase may be locally stabilized by charged impurities in 
the PM. 

Our star ting point is the two-orbital Hamiltonian H = 
Hq + Hi nt 2y, |27|, [28|, l29( with the non-interacting part 
(here written in momentum space) 

H = 2 , [gfc,l^fc.l.cr rf fc,l,g +gfc,2rffc, 2 . CT ^fc,2,g 

fe,<7 

+ £fc,12 (4,1,<A,2,<t + d l,2,a d k,l,<r) (1) 

£k,i = —2t 1 cos(k x ) — 2t 2 cos(ky) — At 3 cos(k x )cos(k y ) 
Sk,2 = —2t 2 cos(k x ) — 2ti cos(ky) — 4t 3 cos(k x ) cos(ky) 
£fc,i = -4t 4 sin(A; a; )sin(A: y ) 

which describes the d xz — and d yz — bands (hereafter "ar" 
and "y") on a square lattice and the interaction part 
includes an intra- and interorbital repulsion and a Hund 
coupling term (here defined in real space): 
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with riiac the occupation operator for orbital a with spin 
a, riia = riiai +n.i a i and Si a the spin operator for orbital 
o. 

In the following we adopt the hopping parameters from 
Ref. [29| and we define the intraorbital repulsion by the 
standard relation U' = U — 3J/2 and J = 0.25C7 for 
definiteness. 

Phases were identified using a fully unrestricted HF 
approximation in large clusters (typically 14 x 14) with 
periodic boundary conditions. Subsequently the energy 
of each solution has been computed in much larger sys- 
tems treated in momentum space. In the latter case we 
neglected a small canting of the MS solution which has a 
negligible effect on the energy. 



In the upper panel of Fig. Q] we show the resulting 
phase diagram of the two-orbital model in the [/-density 
plane at zero temperature. The solid circle is a tricritical 
point: to the left of it the transition is weakly first order 
whereas to the right it is second order. Between the MS 
and the PM we find the new OM phase in which magnetic 
moments at nearest-neighbor sites are at right angles as 
shown in the inset of the upper panel of Fig. [TJ It can be 
seen as a set of toroidal moments [30] mutually antiparal- 
lcl thus it can also be termed an "antiferrotoroidal" state. 
It can also be viewed as the superposition of two magnetic 
stripes perpendicular to each other. The magnetization 
at position R; reads: 



m ( = y^M t exp(Q t .R;). 



(2) 



with Qi = (0, 7r), Q 2 = (n, 0) and Mi, M 2 mutually per- 
pendicular and equal in modulus. The MS is recovered 
by setting one of the M, in Eq. @ to zero. 

The first order character of the transition close to half- 
filling is not surprising given that the model is quite frus- 
trated from the magnetic point of view 3l| . Increasing 
U the system can avoid frustration by keeping in the 
PM phase until the point in which the energy penalty of 
not forming magnetic moments becomes too large driving 
the system through the abrupt transition to a magnetic 
phase. 

We find the topology of the phase diagram shown in 
Fig. Q] to be rather robust although for some parame- 
ters the tricritical point becomes a triple point. Even 
more, we find essentially the same topology in a standard 
one-band Hubbard model supplemented with an explicit 
antiferromagnetic next-nearest-neighbor magnetic inter- 
action J' across the diagonals. The latter frustrates the 
usual Neel state and favors the magnetic stability at Q^. 

The larger stability of the OM phase compared to the 
MS at half filling is due to the fact that this two dimen- 
sional texture opens gaps in the electronic structure in 
both x and y direction whereas the MS in weak coupling 
leaves the direction perpendicular to the stripes metallic. 
Therefore a more stable fully gaped (insulating) solution 
is found for a weaker coupling in the OM case than in the 
MS case. This is illustrated in the lower panels of Fig. Q] 
where we show the density of states (DOS) for both so- 
lutions and U = 0.8eV. For smaller U both solutions are 
metallic but the OM state has deeper pseudogap. 

The dips in the OM DOS close to 1/4 filling and 3/4 
filling can be understood in strong coupling as due to the 
fact that as an electron from the lower (upper) Hubbard 
band circulates a plaquette, the HF potential forces its 
spin to be parallel (antiparallel) to the local direction of 
the magnetization. Thus after a complete loop the elec- 
tron spin gets back to the original direction but the single 
particle wave function acquires a phase e™ = — 1 remi- 
niscent of the staggered flux phase in cuprates[13l |32| |. 
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FIG. 1: (Color online) Upper panel: T — phase diagram of 
the two-orbital model in the [/-density plane. The solid circle 
is the tricritical point. The insets shows the magnetic struc- 
ture of the ordered phases considered in this study. Lower 
Panels: The local HF DOS for the MS and OM states. 



This produces a cone-like dispersion and a semimetallic 
DOS which evolves into the dips in weak coupling. 

More insight on the phase diagram can be obtained 
with the following Landau toy model: 



*/ = sEx-^q) 



(3) 
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It consists of a Gaussian part and a local interaction part 
which we keep up to sixth order to be able to describe 
both second- and first-order phase transitions. We as- 
sume that the susceptibility \ (l) diverges, for a set of 
symmetry related momenta Qj, at an instability line con- 
trolled by the parameter a = x _1 (Qj) and take 71 > 
for stability. The non Gaussian vertexes have been taken 
as momentum independent in the same spirit of Ref. [33j . 
This will be partially relaxed below. 

For /3i > the above model describes an ordinary 
second-order magnetic phase transition as a crosses the 
Gaussian line. For (3i < there is a first-order phase 
transition with a tricritical point at /?i = 0. This can be 
treated analogously to the liquid-solid transition 34 |35| 
and frustrated phase separation[36j. As there, we restrict 
the sums to the Q^ set that in the present case has only 
two elements. Hence we have to choose subsets which 
satisfy the Kronecker delta of the quartic term modulo 



a reciprocal lattice vector. Only two choices are possi- 
ble: either all the (\i are the same and equal to one Q^ 
which leads to the MS or the q^ are equal in pairs to 
one of the Q^. In the latter case we still have to choose 
the relative angle between the two vectorial Fourier com- 
ponents Mi at Qi [c.f. Eq. @] and their magnitudes, 
Mi 37| . For two dimensional textures we find that the 
energy is locally minimized when the magnitude of the 
two amplitudes is the same and the angle between the 
two Fourier components is either tt/2 (which leads to the 
OM state) or zero. The latter is a new phase in which 
the real space magnetization is zero in one sublattice and 
forms an AF structure in the other sublattice. This is the 
only phase we have found in which the modulus of the 
magnetization is not uniform. Since the charge density is 
a scalar, it couples with the square of the magnetization. 
It follows that this phase will have both spin and charge 
order (SCO) reminiscent of the charged stripe phases in 
cuprates[8|,lll(- In the present case the charge ordering 
wave vector is (jr, 7r). It is suggestive that this phase is 
similar to the charge ordering in BaBi03 which becomes 
superconducting when doped[6|, |7j. In the following, for 
simplicity, we neglect the effect of the charge relaxation 
upon the energy which, in any case does not affect the 
quadratic terms. Indeed we find in HF that this effect is 
small. 

Calling Mt the magnitude of the magnetization in real 
space, we find that the MS and the OM phase have the 
same energy [38j| given by SJoalms — aJW|/2 -I- j3\M^, + 
7iMj. while for the SCO phase Sfsco — %m,ms/2- 
These energies are easily rationalized from the fact that 
the non Gaussian terms in Eq. (|3|) are local in real space 
and in the SCO state only half of the sites are magnetized 
(c.f. inset in Fig. [J). 

The phase diagram is shown in Fig.^a). The dashed 
line is the limit of metastability of the magnetic phases. 
The insets show the behavior of / in the metastable re- 
gions. The 1/2 factor in the energy makes the MS and 
OM configurations more stable than the SCO configura- 
tions in the ordered region. As soon as the phase bound- 
ary is crossed and one reaches the PM phase the situation 
is reversed. The lowest energy phase above the PM is the 
"hidden" SCO phase. 

In the restricted HF approximation for the param- 
eters of Ref. [29( and positive doping we have found 
that the PM-SCO is second order, thus it can not be- 
come metastable as in the Landau toy model but, of 
course this is sensitive to the parameter set. For other 
parameters 28j we have even found a region of SCO phase 
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in the phase diagram at high doping. 

The Landau toy model phase diagram is of course over- 
simplified and does not even reproduce the full physics 
of the HF phase diagram. One can construct a canoni- 
cal Landau theory as follows. First the magnetization is 
written as the sum of products of slowly varying parts 
Mj(R) times rapidly varying parts exp(Qi.Rj). The en- 
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FIG. 2: (Color online) Phase diagram for the Landau toy 
model (upper panel) and the full Landau model (lower panel) . 
The dashed line is the limit of metastability of the ordered 
phases. The limit of metastability of the PM phase is given 
by the line a — 0. The thick lines represent weakly first 
order transitions ending at the tricritical point. The insets 
in the upper panel show the behavior of the energy for the 
magnetic phases (full line) and the SCO phase (dashed line) 
in the regions of metastability. 



ergy can be expanded close to the Gaussian line in terms 
of the invariants generated by the Mj(R) and its gra- 
dients. Here we are interested in uniform phases so we 
neglect the spacial dependence of the order parameters 
Mj reaching Eq. ©. 

We have two second-order invariants Ii = Mj.Mj. The 
fourth-order invariants are powers of the second-order 
ones and I 3 = (Mi.M 2 ) 2 and I A = (Mi x M 2 ).(Mi x 
M2). The latter can be eliminated since it depends on 
the other invariants through the relation 1^+1^ = ii/ 2 . 



The energy can be written as: 



SF = -a^M l .M J + /3 1 ^(M. i .M l ) 2 +/3 2 (M 1 .M 2 ) 2 

i=l i=l 

+ /3 3 (M 1 .M 1 )(M 2 .M 2 ) 

2 2 

+ 7l ^(M 4 .M l ) 3 + 72 (M 1 .M 2 ) 2 ^M I .M l (4) 

i=i i=i 

+ 73 [(M 1 .M 1 )(M 2 .M 2 ) 2 + (M 1 .M 1 ) 2 (M 2 .M 2 )] 



In this case the energy of the different phases is: 
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with Si = (3 U B 2 = (2/?i+/? 3 )/4, £ 3 = (2/3i+/3 2 + /3 3 )/8, 
d = 71, G 2 = (71 + 73) /4, G 3 = (71 + 72 + 73)/16. 

The parameters of the model can be fixed by compar- 
ing observables with experiment. One can consider the 
former to be function of two or more control parame- 
ters (like U, n, P and T). Taking two control parame- 
ters for simplicity, each phase considered alone can have 
a different tricritical point with the paramagnet. If in- 
stead all phases are allowed to compete, depending on 
parameters one can have more complex phase diagrams 
including triple points. Given the local character of the 
electronic interactions, we expect the parameters to be 
not far from those of the Landau toy model Eq. © , which 
corresponds to flj = /3i , G h = 71 . Indeed one can obtain 
the same topology as in HF taking B\ — /3i, B 2 — 1.07/?i, 
B 3 > max(0,/3i/2), G 2 = I.I71 and G x = G 3 = 71 [c.f. 
Fig. [2](b)]. 

The four phases identified (PM, OM, MS, SCO) are 
the only phases allowed by the model Eq. ([4]) thus our 
study is exhaustive close to the Gaussian line (i.e. close 
to the tricritical point in the first order region). The 
exception are special combinations of parameters where 
a family of solutions become degenerate as in the Landau 
toy model [38|. Canting of the magnetic moments respect 
to the identified solutions require higher order terms in 
the energy. 

It is easy to check that the structure factor averaged to 
take into account a superposition of randomly oriented 
domains is the same for the OM and the MS phases, thus 
in principle it can be difficult to distinguish among them 
with a magnetic neutron scattering experiment alone on 
polycrystalline samples. In practice, nuclear neutron 
scattering detects also a lattice distortion which breaks 
the G4 symmetry of the lattice and probably stabilizes 
the MSJ24I31J. It is possible that the OM state can be 
stabilized suppressing the structural transition by pres- 
sure or chemical substitution. 

The SCO phase may be locally stabilized by charged 
impurities in the paramagnet. We expect the concomi- 
tant spin and charge order to appear around charged 
nonmagnetic impurities which may be observed with lo- 
cal probes like NMR, NQR and scanning tunneling spec- 
troscopy. 

In conclusion we have exhaustively analyzed possible 
magnetic phases competing with superconductivity in 
FeAs layers close to the Gaussian instability line of the 
paramagnet. Because of magnetic frustration we find a 
tendency for PM-magnetic transitions to be first order 



as shown in the HF approximation. This makes non 
Gaussian terms relevant and generates an anharmonic 
order parameter as shown in HF and using Landau the- 
ory. Contrary to general beliefs the two-orbital minimal 
model does not show the MS phase close to half-filling. 
Instead we find a new phase in which magnetic moments 
acquire an unusual orthogonal configuration. We find 
another low energy phase with spin order and charge or- 
der at momentum (n, n) which provides an obvious link 
among charge fluctuations, possibly relevant for super- 
conductivity, and magnetism. Our result provide a guide 
of likely configurations to be found in the phase diagram 
of layered FeAs based compounds and a Landau frame- 
work to study them. 
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Notice that because of spin rotational symmetry the en- 
ergy is invariant with respect to a global rotation of the 
magnetizations. 

For the Landau toy model we find that the energy of 
solutions with Mi and M2 orthogonal depend only on 
Mt = Mi + Mi and not on 5M 2 = M? - Ml leading to 
a family of degenerate phases. 



